Screening and inplane magnetoresistance of anisotropic two-dimensional gas 
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In order to split the influence of the orbital and spin effects on the inplane magnetoresistance of 
a quasi two-dimensional (2D) gas we derive its linear response function and dielectric function for 
the case of anisotropic effective mass. This result is used for the calculation of elastic transport 
relaxation time of a quasi two dimensional system in a parallel magnetic field. The relaxation time 
is proved to be isotropic in the low density limit for the case of charged impurity scattering, allowing 
to separate the two contributions. 
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Among a variety of experimentally used two di- 
mensional semiconductor structures some possess the 
anisotropy of Fermi surface originating from that of a 
bulk material, including n-type Si-MOSFET's on other 
than (100) surfaces 0, AlAs heterostructures Q, p-type 
Si and GaAs structures 0- The transport properties 
of such anisotropic semiconductors has been addressed 
largely |J. LJ. L| an d it is well known, that the Fermi sur- 
face anisotropy gives rise in general to the anisotropy 
of relaxation time, even if the scattering potential is 
isotropic 0. 

Another possibility is the externally introduced 
anisotropy through the application of parallel to a quasi 
2D layer magnetic field, known to deform a Fermi sur- 
face due to the so called orbital effect @ . In this case, 
however, there is additionally a coupling of parallel field 
to the particles' spins leading to the partial spin polar- 
ization of the system [tj [n| and subsequent change of 
screening flal (so called spin effect}. Although substan- 
tial interestH H OH 13 EH 111 has been recently 
paid to the longitudinal magnetoresistance (MR) studies 
of quasi two dimensional systems, it remained unclear so 
far how the two effects play together [Toj . 

In this letter we address the screening properties of 
the anisotropic 2D gas deriving the linear response func- 
tion for the case of elliptic Fermi surface. This result is 
then used to evaluate the transport relaxation time r for 
elastic charged impurities scattering, which, surprisingly, 
turns out to be isotropic in the low density limit. Further- 
more, for partially spin polarized anisotropic 2D system 
the relaxation time for each spin subband is also shown 
to be isotropic, allowing finally to reach the main result 
of the paper — to separate the influence of the orbital 
and spin effects on the longitudinal magnetoresistance of 
a diluted quasi 2D gas. 

In the following we utilize the simplest possible rep- 
resentation of anisotropy — the elliptic Fermi surface in 
the effective mass approximation. The interactions are 
treated in the random phase approximation (RPA), so 
that the screened linear response function equals to lin- 
ear response of a free particles gas [lq]: 
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where E^ is kinetic energy of a quasiparticle with mo- 
mentum k, /o(-E'k) - the zero temperature Fermi-Dirac 
distribution function, g v ,g s — the valley and spin de- 
generacies. The excited state \n) contains a single pair 
of a quasiparticle with momentum k + q and a quasihole 
with momentum k, and the matrix element of a den- 
sity fluctuation operator (pq) is equal to unity. The 
last property originates entirely from the Bloch type of 
Hamiltonian eigenfunctions in the effective mass approx- 
imation, similar to the isotropic case |l6j| . 

We perform the following change of coordinates to 
rewrite the integral in a spherically symmetric form: 
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The kinetic energy depends solely upon the length of 
the distorted wave vector k — E^ = Ep((k x /k x ) 2 + 
(ky/ky) 2 ) = Ep |k| 2 and for integral Q we have: 
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In view of spherical symmetry the integral value in @ 
is invariant to rotation of vector q, hence the linear re- 
sponse function depends upon its length q only. The 
integral in (0 would be the same for isotropic Fermi sur- 
face as well, thus the only difference from Stern's linear 
response function [l7| is due to the normalizing prefactor 
proportional to the density of states of anisotropic Fermi 
gas D = g v g s ^m x m y /2'Kh 2 . Finally, we get for \'- 
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The linear response function |0J of anisotropic system 
depends upon the direction of perturbation wave vector 
q through the variable q = {{q x /k x ) 2 + {qy/ky) 2 ) 1 / 2 , so 
that the screening does become anisotropic, in contrast 
to the case of isotropic Fermi surface 01 • Note, however, 
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that this anisotropy is the same as that of kinetic energy 
as a function of momentum, since q 2 = E^/Ep. 

At the same time we are able to find the RPA static 
dielectric function |l6j : 



s(q,Q) = l-V(q)- X SCr (ci,Q) = l 
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where V(q) = 2ire 2 /q is the 2D Fourier transform of the 
bare Coulomb interaction potential, Qtf = 2ire 2 D — the 
Thomas- Fermi screening parameter 

We now turn to the application of our results to the 
transport properties of anisotropic 2D Fermi gases. Ex- 
pression (J5J for dielectric function enables one to find 
the elastic transport relaxation time. In general r(k) is 

anisotropic and calculation of this could be a rather com- 
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plicated procedure 6, lq. As we show below, for the case 
of screened charged impurity scattering, t is still isotropic 
in the low density regime and can be obtained analyti- 
cally. One has for elastic scattering transport relaxation 
time: 
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where v^. = h 1 dE-^/d\i is the particle's group veloc- 
ity. Note that the relaxation time isotropy is already 
implicit in this expression and is verified in what fol- 
lows. First term in the integrand is the square of a 
scattering matrix element in the Born approximation Q , 
where we have neglected all the form-factors of the 
real Coulomb interaction between a quasi 2D electron 
(hole) and a charged impurity 0. Applying the change 
of coordinates J2J to the k'-space we find from (JSJ: 
V(g)/e(q, 0) = 2ire 2 (q + qTF^iQ))^ 1 which in the low 
density limit (k^ ,ky <§; qrF) reduces to pj|: 

V(q)/e(q,Q) = 2ire 2 q T 1 F( f>- 1 (q) (7) 

The second term in © accounts for a loss of initial 
velocity in a scattering event and is similar to the factor 
1 — cos(k, k') in the isotropic case Changing the 

coordinates to polar ones k x ,k y — » k, 9 we rewrite this in 
terms of a scattering angle 9: 1 — Vk ■ /v 2 . = 1 — cos 9 + 
A(k) sin 9 , where number A(k) depends upon the initial 
particle momentum k. Hence in the low density limit we 
find: 
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since q — (2 — 2cos6') 1 ^ 2 , and the odd in 9 part of the 
integrand gives no contribution to the integral. This fi- 
nal expression for the elastic scattering time in the low 
density limit is essentially the same as in the isotropic 
case, the only difference represented by the reciprocal 
density of states D -1 in the prefactor. This finding has 
two important consequences: in the low density limit, 
within the elliptic deformation of the Fermi surface, the 



transport relaxation time (i) remains isotropic andfii) 
increases proportionally to the density of states D [2(| • 

Comparing to previous studies we find that our result 
for the relaxation time isotropy recovers the one derived 
earlier for the short range scatterers 6], since for the 
case of zero spin polarization considered so far the rele- 
vant Fourier components of the screened impurity poten- 
tial J7J do not depend on wave vector at all, according 
to The predicted increase of the relaxation time 

with the Fermi surface deformation implies of course the 
increase of conductivity, i.e. the negative magnetoresis- 
tance in the parallel to field direction caused by the or- 
bital effect, as the effective mass in this direction remains 
unchanged 0. This is in contrast to positive MR found 
for orbital effect in Ref. , where the change of screening 
has been neglected. Our prediction could be easily ver- 
ified in the experiments on wide quantum wells, where 
the inplane magnetoresistance is mostly due to the or- 
bital effect. 

As was mentioned earlier, apart from the Fermi sur- 
face deformation the parallel magnetic field couples to 
particles' spins, resulting in a partial spin polarization 
of a system. Similar to the isotropic case [ijj there is 
no more a single Fermi surface, but different ones for 
different spin projections onto magnetic field axis. We 
would like to treat only the spin conserving processes, 
which means that the linear response of such partially 
polarized system is simply a sum of the responses from 
different spin subbands 15]. Similar to the case of un- 
polarized system (expression (@J), the anisotropy of the 
linear response of partially polarized system is again the 
same as that of the kinetic energy as a function of wave 
vector q, which leads to the isotropy of the relaxation 
time for each spin subband in the limit of low density, as 
we show below. 

Calculating the transport relaxation times one 
should write the integrals of type 10 for major (f ) and 
minor (|) spin subbands separately. In the low den- 
sity limit, when the Fermi wave vectors of both sub- 
bands satisfy k%''^, ky''^ <C qrF, the Fourier image 
of the screened impurity potential V(q)/e(q) has the 
anisotropy of the linear response function and kinetic en- 
ergy, according to Q. This means that the change of co- 
ordinates @ applied to the integrands leads to the same 
expressions for as one gets in the isotropic case |l5j . 
Thus the relaxation times are isotropic, the only effect 
of anisotropy been again to normalize the absolute value 
of t*" through the density of states dependent prefac- 
tor (JSJ) . We are now able to write down the conductivity 
tensor of the partially polarized anisotropic system: 
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where a\ are the major and minor spin subbands con- 
ductivity tensors, e, 715, To, Dq — respectively, electron 
charge, density of 2D particles and the zero field isotropic 
relaxation time and the density of states. D and m are 
respectively the density of states and the effective mass 
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tensor in magnetic field. The last term in © stands for 
the Dolgopolov-Gold's calculated change of the conduc- 
tivity of isotropic system as a function of its degree of spin 
polarization £ = (n-f — n.j,)/(nj + rii) 15]. Note that for 
the case of a quasi 2D system with zero field anisotopic 
mass this result is valid only for a parallel field applied 
along the main axes of symmetry, otherwise the Fermi 
surface loses its ellipticity in magnetic field. Let us show 
finally how this simple expression allows to separate im- 
mediately the contributions from spin and orbital effects 
on the longitudinal magnetoresistance of a diluted quasi 
2D system. 

We focus on the recent magnetoresistance studies 
of the 2D electron gas of AlGaAs/GaAs heterostruc- 
ture [1(1 ]. This system is isotropic in zero magnetic field, 
hence, according to © the anisotropy of experimental 
MR [Hj is due to the effective mass change in a perpen- 
dicular to magnetic field direction 

mx/mo = px(B)/p\\(B) , 

where _L and || mark the resistances measured in perpen- 
dicular and parallel to inplane magnetic field directions, 
respectively. The spin effect contribution in J^J is given 
thus by: 



1 p±(B)- Pll (B) 
px(B = 0)-p||(B = 0) 
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The degree of spin polarization depends upon both the 
Zeeman energy gp,sB and the effective mass at a given 
field value: £(B) = {E F] - E Fl )/{E F1 + E Fl ) = gp, B B ■ 
D/ns, where Epf y are the kinetic parts of the Fermi en- 
ergy for two spin subbands, g, p, B — the Lande factor and 
Bohr magneton. Equivalently the last equation reads: 
i{B) = B/B p ■ y/m ± (B)/mo, where B p = 2E%/gp B is 
the full spin polarization field in the absence of orbital 
effect |T|. 

In the inset to Fig^we show the effective mass growth 
extracted in the above manner from the p±,p\\ data of 
Ref. ^3 • The effective mass in perpendicular to field di- 
rection grows by about 30% in moderate fields, as caused 
by the orbital effect Q . For such a slight deformation the 
utilized approximation of elliptic Fermi surface should 
work reasonably, in contrast to general case Q . The par- 
tial spin polarization has major effect on MR, as reflected 



by its small anisotropy ^} , and leads to a roughly three- 
fold resistance increase as is shown in the body of Fig^ 
The full spin polarization has not been reached in [l£| 
that's why the saturation |15| of the geometrical mean of 
the parallel and perpendicular to inplane field resistivities 
predicted by fTH| is not seen in Fig^ Fitting to calcu- 
lation 0] we obtain the B P value 17.3 T which corre- 
sponds to a Lande factor g* w 2.1 at ns ~ 3 x 10 11 cm~ 2 , 
in agreement with previous studies |lCl ll 1| . 



We would like to add a note here, concerning the ap- 
plicability of the MR data analysis presented above to 
the real interacting quasi 2D systems. Apart from the so 
far considered single particle effect of Fermi surface de- 
formation, the squeezing of the 2D layer by the parallel 
magnetic field additionally changes the form-factors of 
Coulomb interaction between particles [I| and increases 
the Wigner-Seitz ratio [y|, which can in principle lead 
to the renormalization of the zero field effective mass and 
g- factor [21|. In presence of such many-body effects the 
spin effect contribution cannot be extracted with for- 
mula H10[). Experimentally this means that the geomet- 
rical mean of resistivities yj p±(B) ■ p\\{B) does not satu- 
rate upon the reach of full spin polarization. The inplane 
magnetoresistance anisotropy, however, should still give 
the anisotropy of effective mass pj_/p\\ = m±/mu, similar 
to the single-particle picture, although the independent 
measurement is required to find the full spin polarization 
field B% [30. 

In conclusion, we have derived the linear response 
function and dielectric function of 2D Fermi gas with 
anisotropic effective mass. In the low density limit the 
screened charged impurity potential is shown to possess 
the same symmetry as the kinetic energy as a function of 
wave vector. As a result the elastic transport relaxation 
time t is isotropic in this limit, even if the 2D system is 
partially spin polarized. This finding allows us to sepa- 
rate for the first time the influence of the orbital and 
spin |l5j ] effects on the inplane magnetoresistance of a 
diluted quasi two dimensional system. 
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FIG. 1: In the inset the effective mass growth caused by the 
orbital effect [81 is shown, extracted from the longitudinal MR 
data of Ref. [lfJl (ns ~ 3 x 10 10 cm -2 ) as described in the text. 
The spin effect contribution is shown in the figure body. The 
fit to theoretical dependence [l5| (dashed line) gives full spin 
polarization field in the absence of orbital effect B% « 17.3T 
corresponding to a Lande factor g* ~ 2.1 



